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Abstract 

Wc compute dilcpton invariant mass spectra from the decays of p mesons produced by photon reactions 
off nuclei. Our calculations employ a realistic model for the p photoproduction amplitude on the nuclcon 
which provides fair agreement with measured cross sections. Medium effects are implemented via an earlier 
constructed p propagator based on hadronic many-body theory. At incoming photon energies of 1.5-3 GeV 
as used by the CLAS experiment at JLAB, the average density probed for iron targets is estimated at about 
half saturation density. At the pertinent 3-momenta the predicted medium effects on the p propagator are 
rather moderate. The resulting dilepton spectra approximately agree with recent CLAS data. 
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1. Introduction 

The investigation of hadron properties in hot 
and/or dense matter is of fundamental interest in 
the context of approaching the transition(s) into 
a chirally restored and/or deconfined plasma of 
quarks and gluons. Intriguing effects have been ob- 
served in dilepton spectra measured in high-energy 
heavy- ion collisions pj, t2] , which are consistent with 
a strong broadening of the p-mcson spectral func- 
tion by about a factor of ^^3 in hot and dense 
hadronic matter Q • The medium modifications 
of the p are believed to be largely driven by the 
baryonic component of the medium. The rapid ex- 
pansion of the fireball formed in heavy-ion reac- 
tions implies that the emission spectra encode a 
rather large range of temperatures and densities 
of the evolving medium. It is therefore desirable 
to test the medium effects in a static environment, 
such as provided by ground-state nuclei. Hadronic 
models predict appreciable medium effects in cold 
nuclear matter, e.g., an increase of the width of 
low-momentum p mesons at saturation density by 
a factor of 2-3 over its vacuum value Isl. Several 
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experiments have recently been conducted to mea- 
sure p production off nuclei, in both proton- Q and 
photon-induced [1, 0, S| reactions. In Ref. 0] a 
rather small dilepton signal for the p has been re- 
ported with a mass distribution compatible with a 
dropping mass, while Refs. @ and 0] extracted a 
moderate broadening with little, if any, mass shift 
in 7r"'"7r~ and e"'"e~ mass spectra, respectively. An 
inherent feature of nuclear production experiments 
is that rather large projectile energies are required 
to supply the rest mass of the p. These impart an 
appreciable 3-momentum on the p meson relative 
to the nucleus which enhances the probability for 
decays outside the nucleus thus reducing the effec- 
tive density probed by these experiments. Never- 
theless, valuable constraints for cold nuclear matter 
effects and their 3-momentum dependence on exist- 
ing models for in-medium p spectral functions can 
be expected. 

As in heavy-ion reactions, dileptons are of spe- 
cial interest due to their negligible final-state in- 
teractions. However, the initial states in heavy- 
ion collisions and nuclear production experiments 
are quite different. In the former case, the sim- 
plifying assumption of a thermal heat bath can be 
made, while in the latter case a reliable description 
of the elementary production process is mandatory. 
The pertinent baseline reaction on a single nucleon, 
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7 -/V — > e+ e N , has been studied in several theo- 
retical works 0, [13, 



[13, El- Generally, the low- 



energy cross section is dominated by baryon reso- 
nance formation, while at photon energies of ~1.5- 
2 GeV f-channel exchange processes are expected 
to take over. Applications to nuclear targets can 
be found in Ref. where a schematic model for 
the production process has been implemented into 
a transport simulation for final-state interactions, 
as well as in Ref. for the CLAS data 0,|3l. 

In the present work we combine a microscopic 
model for p photoproduction on the nuclcon [12l | 
with an in-medium p spectral function computed 
in hadronic many-body theory [3, [lE] ■ The pro- 
duction model is largely based on meson/Pomeron 
exchange which properly accounts for the cross sec- 
tion above photon energies of 2 GeV. At lower en- 
ergies, we supplement additional s-channel reso- 
nance excitations with parameters directly taken 
from the in-medium selfenergy of the p spectral 
function [3, [l^. This establishes consistency be- 
tween the production process and in-medium ef- 
fects, and leaves no additional free parameters for 
the resulting cross sections for p production and 
dilepton invariant-mass spectra. 

In Sec.[2]the p photoproduction model on the nu- 
cleon is presented and checked against total cross 
sections and dilepton invariant-mass spectra for 
deuteron targets. In Sec. Owe apply our model to 
dilepton spectra off nuclei utilizing the in-medium p 
spectral function at densities estimated from the de- 
cay kinematics corresponding to the incoming pho- 
ton spectrum in the CLAS experiment. We finish 
with conclusions in Sec. [H 

2. p Photoproduction on the Nucleon 

Our starting point is the photoproduction arnpli- 
tude for 7p — > e"*" e~ p developed by two of us [12j . 
It accounts for cr, /2, 2tt and Pomeron i-channel 
exchange as well as nucleon s- and it-channel pole 
contributions, and gives a good description of ex- 
perimental cross sections at photon energies above 
~2 GeV, cf. the dotted line in Fig. [TJ (A similar 
model in Ref. employs a stronger cr-exchange. 



see Ref. [12|,[16| for further comparison.) At smaller 
photon energies, as part of the photon beam used 
by CLAS at JLAB (go ^ 1-3.5 GeV), baryon reso- 
nances are expected to become important 17| . Here 
we adopt the same set of resonances as used in 
Ref. [ij] to describe total photoabsorption spectra 
(to constrain the in-medium p spectral function) 
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Table 1: Resonance parameters (columns 2, 5 and 6) used in 
the elementary photoproduction amplitude. The dimension- 
less coupling constants (all other quantities are in [MeV]) 
are fixed to reproduce the values for the total and partial 
vacuum on-shell widths (columns 3 and 4) as in Refs. |15|| 
(the cut-ofi' values of Refs. [3I. Il5l| are kept fixed). 



summarized in Tab. [10 Due to the rather large 
photon energies involved we employ relativistic in- 
teraction vertices defined by the following p-N-B 
Lagrangians: 
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where (t^'' = f[7'^ ,7l, = ft" " 

Isospin, spin and parity of the resonance B are de- 
noted by IJ^ , and Tj, Ti are the usual isospin-1/2, 
1/2 3/2 transition matrices[l As in Ref. 14|. 
we utilize an improved version of the vector dom- 
inance model (VDM) [3, [13] which allows for a 
direct ^-N-B coupling and thus a better simulta- 
neous description of hadronic and radiative decay 
widths (the j-N-B coupling follows by replacin g Pn 
with Af^ in Eqs. HJ; the parameter = 0.7 ]l3| 
controls the deviation from naive VDM). The pa- 
rameters (listed in Tab. [T]) are adjusted to recover 
the same partial decay width for B pN as the 



^Thc pNN formfactor has been reduced to 600 MeV im- 
proving consistency with Ref. [T3l . together with a 10% re- 
duction of the p-N-N{1720) coupling constant. The extra 
coupling to the A'^(2090) increases the p production cross 
section around photon energies of ~2 GeV by ~15%. 

^The spin-5/2 resonance is treated in a simplified way as 



Ref. [l|] 



state with amended spin factor. 



(updated) values in Rcf. [Sj (based on Ref. [ij]), 
including monopole formfactors 



A2 



(2) 



with cutoff parameters ApSN The formfactors 
are consistently evaluated in the laboratory frame 
with q the three-momentum of the incoming photon 
(or p in the nuclear rest frame). The spin- 1/2 and 
-3/2 baryon propagators are, respectively, taken as 
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with masses and total widths Tb as given in 
Tab. [TJ We furthermore assume a linear increase of 
the in-medium resonance widths with density as in 
Refs. 0, We have verified that neglecting the 
momentum dependence of in the propagators 
has an insignificant impact on our results. Includ- 
ing s- and M-channel graphs, the baryon resonance 
parts of the p-production amplitude take the form 
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for intermediate spin-1/2 and -3/2 states, respec- 
tively; p and p' are the momenta of the in- and out- 
going nucleon, q (k) is the photon (p) momentum, 
and X/ = 2 (4/3) an isospin factor for 7=1/2 (3/2) 
resonances. The vertices T follow from Eq. ([T]): 
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for negative parity resonances. 

It is now straightforward to implement the pro- 
duction amplitude of Ref. [l2|, augmented by 
baryon resonances, into a mass differential cross 
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Figure 1: Total cross section for 7p — » p p" 
of incident photon energy in the laboratory frame. The 
calculations are based on either a fixed p-mass of 770 MeV 
(dash-dotted and dashed line) or a full vacuum spectral func- 
tion (dotted and solid line), either with (dashed and solid 
line) or without (dash-dotted and dotted line) baryon reso- 
nances. Data are from Refs. |23L |24||. 



section per nucleon for exclusive p and dilepton pro- 
duction. For the latter one obtains 
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where follows from summing A^f over the pho- 
ton polarization in the elementary processes, 

T'^{q,p,k)^ J2 Mriq,P,k)e,{q), (8) 

ie{a,4,,N,h,B} 

and the dilepton final state is represented by 
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Eq. ([7|) contains an average over the Fermi motion 
of the incoming nucleon (pg = p ^ + m^) as needed 
for nuclear targets in the next section. Further- 

L IT 

more, Gp denote the longitudinal and transverse 
components of the electromagnetic correlator [3] 
(in naive VDM, one has Gp^^ = [mf^ f j g^D^p''^ 
where D^'^ is the p propagator), and P^^;^ are pro- 
jection operators. 
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Figure 2: Differential p-meson pliotoproduction cross sec- 
tion. Calculations with (solid lines) and without (dashed 
lines) resonance contributions are compared to data |23|I . 
with 2 curves each representing the upper and lower end 
of the experimental photon energy window. In each case the 
full vacuum p-meson spectral function has been used. 



Our main interest in the present paper concerns 
the shape changes in the dilepton mass spectra in- 
duced by the in-medium p propagator encoded in 
Gp^^ . Note, however, that Eq. ([7]) also accounts 
for the reduction in dilepton emission due to ab- 
sorption of the p meson propagating in a nuclear 
medium at fixed density, via the in-medium reduc- 
tion of |Gp^^p. For finite nuclei, this effect causes 
an appreciable decrease of the total e+e~ produc- 
tion cross section, relative to a simple scaling with 
nuclear mass number, A (referred to as nuclear 
transparency ratio, Ta] see, e.g., Refs. [20l. [ill. [2^ 
for the cases of u and </> pliotoproduction). 

We first test our production amplitude in the pro- 
cess 7p — *■ p'^p. The photon-energy dependence of 
the p-production cross section is shown in Fig. [TJ 
The contribution of the resonances nicely fills in 
low-energy strength that was missing in the orig- 
inal model of Ref. and becomes negligible at 



Figure 3: Dilepton invariant-mass spectrum for p photopro- 
duction off deuterium with (solid line) and without (dashed 
line) baryon resonance contributions, compared to CLAS 
data after subtraction of u! and contributions 0, H ■ 



energies beyond 2.5 GeV. Our calculations also il- 
lustrate that the inclusion of the free p width (as 
given by the vacuum spectral function of Ref. [lj|) 
further improves the agreement with the low-energy 
cross section (however, in this region the extraction 
of the data is beset with significant model depen- 
dence (Hi)- The scattering-angle differential cross 
section (Fig. [2|) reveals that the resonance excita- 
tions provide large contributions at large scattering 
angle which is supported by experiment. 

Next we apply our model to dilepton invariant- 
mass spectra off deuterium. To mimic a finite 
nucleon-momentum distribution and rescattering 
effects in Gp we use a small average density of 0.1 giQ 
in Eq. ([7]) (folding over a realistic density distri- 
bution gives similar results). The incoming pho- 
ton energies are weighted in 6 bins from qq = 1- 
3.5 GeV to simulate the Bremsstrahlungs-spectrum 
used by CLAS [2^. The shape of the e+e" spec- 
tra [1] is reasonably well reproduced, except for 
masses above 0.9 GeV where additional production 
processes might become relevant, see also Ref. Q. 

3. Dilepton Spectra off Nuclei 

To evaluate medium effects for nuclear targets, 
we first have to estimate the densities probed for a 
given nucleus. If the p meson were produced at rest, 
the density at its creation point would be a good 
approximation. However, since we are considering 
rather high photon energies the p meson will travel 
a significant distance before it decays. Based on the 
assumption that the (medium effect on the) p in- 
stantaneously adjusts to the surrounding medium, 



the relevant density for the dilepton spectrum is the 
local density at the decay point, which we estimate 
as follows. For the incoming photon the interaction 
point is distributed according to a Woods-Saxon 
density profile (weighted by volume). The average 
travel distance of the p from its production to decay 
point is then calculated as 

-1/2 



(10) 



where \v\ is the p three- velocity and t its aver- 
age lifetime (time dilated by a Lorentz 7 factor). 
Under the present conditions, the latter is roughly 
~ 1 fm/c from the underlying in- medium spectral 
function, cf. Fig. [H The velocity is estimated from 
the incoming photon energy for an on-shell p in 
the limiting case of forward production where the 
bulk of the differential cross section is concentrated 
(recall Fig. ^ where nuclear Fermi motion is ne- 
glected). The travel length L obtained in this way 
is then integrated over all production points result- 
ing in the following distribution of decay points at 
a given density g^, 

Nigcc) = fg{r, 0, e) sm{e) 6{g{r, L, 6) - g,) d\, 



g{r,L,9) 



90 



1 -I- exp 



,(11) 



with z = 0.55 fm and c = 4.05 fm for iron. At an av- 
erage incoming photon energy of ~2.1GeV (repre- 
sentative for the CLAS experiment [2^) the average 
density at the decay point amounts to 0.5 go. Vary- 
ing the photon energy between 1.5 GeV and 2.5 GeV 
affects the average density by about ±0.1 go- We 
therefore display our dilepton spectra on iron for 
a density range of gN = 0.4-0.6 £>oIl Note that 
lower photon energies (probing larger densities) im- 
ply a smaller p-meson phase space which is there- 
fore biased toward lower invariant masses. While 
the incoming photon energy spectrum is properly 
included at a given density via Eq. ^ , the density- 
energy correlation is neglected. However, across the 
above range, the density variation of the dilepton 
spectra turns out to be quite moderate. A more 
accurate evaluation of this correlation should also 
include an exponential decay distribution of the p 
decays in L. 



^This also reflects some of the uncertainty introduced by 
letting all p mesons decay after a flxcd distance L, as com- 
pared to a distribution in L, since a different photon energy 
translates into a different L. 




Figure 4: fn-medium p-meson spectral function at various 
3-momenta and nuclear density 0.5 go Il5l : solid (dashed) 
lines: transverse (longitudinal) parts (identical at ij = 0). 



The main in-medium input to Eq. ([7]) is the p 
spectral function of Refs. IJ, ll5[ which is displayed 
in Fig. [3] for transverse and longitudinal modes at 
various 3-momenta and at = 0.5 ^iq- At 3- 
momenta relevant for CLAS {q ~ 1-2 GeV) the 
medium effects are significantly reduced compared 
to ? = (a consequence of the typical formfactor 
cutoffs, ApBN — 0.6 GeV; the reduction is more 
pronounced than, e.g., in the spectral function of 
Ref. [l^ due to larger formfactor cutoffs used in 
there). In addition, a noticeable difference be- 
tween longitudinal and transverse modes develops, 
the latter exhibiting an upward mass shift which 
is due to both pion cloud and P-wave resonance 
excitations. Note that in applications to dilepton 
spectra at CERN-SPS the in-medium spectral func- 
tion is pred ominantly probed at 3-momenta below 
isl . 3l- This reiterates the notion that the 



1 GeV 

CLAS data provide a novel test of the spectral func- 
tion at high 3-momentum. 

Our results for the dilepton invariant-mass spec- 
tra on iron are compared to the CLAS "excess" 
spectra in Fig. [5] using the density range as esti- 
mated above. For each density, the normalization 
is adjusted to the integrated strength of the data. 
Alternatively, one can determine the normalization 
by a least-square fit resulting in x^/^=l-29 (1.4) 
per data point (not) including the resonance con- 
tributions in the production process, compared to 
X^/iV=1.34 (1.49) when normalizing to the data. 
In either case, the agreement with the data is fair 
(the in-medium broadening of the nucleon reso- 
nances has very little impact on the dilepton spec- 
tra). A slight discrepancy with the data for masses 
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Figure 5: Theoretical calculations of dilepton spectra for 
photoproduction off iron compared to CLAS data The 
bands represent the nuclear density range = 0.4-0.6 go- 
The curves in the upper (middle) panel are calculated with- 
out (with) baryon resonances (Tab.^ in the elementary pro- 
duction amplitude, while the lower panel shows the iron-to- 
deuteron ratio for the full calculation. 

of AI = 0.8-0.85 GcV may allow for a small at- 
tractive mass shift of about -15MeV0 Overall, the 
rather moderate medium effects in the (transverse 
and longitudinal parts of the) p spectral function 
at high 3-momentum (as seen in Fig. H]) are es- 



sentially in line with the CLAS spectral^ There 
are further effects which could modify our spec- 
tra at the several percent level, e.g.. in- medium lo- 
meson decays along with interference/mixing with 
the p. a and /2 tadpole diagrams are not included 
in the spectral function; implementing the former 
with the coupling strength employed in the elemen- 
tary production process [l^ generates an attrac- 
tive mass shift of about —10 MeV for the iron tar- 
get. Further processes in inclusive production (e.g., 
'y + N p + N + Tr), additional resonance strength 
to accommodate large angle-scattering at high pho- 
ton energies (recall lower panel in Fig. [5]), or a 
more elaborate treatment of the baryon-resonancc 
widths, might also play a role. 

4. Conclusions 

We have performed an essentially parameter-free 
calculation of p photoproduction off nuclei, combin- 
ing a realistic model for the elementary production 
process with a hadronic many-body spectral func- 
tion [ist which was extensively used before in the 
interpretation of dilepton spectra in heavy-ion colli- 
sions. An earlier constructed photoproduction am- 
plitude has been supplemented with resonance 
contributions as implicit in the in-medium p spec- 
tral function. A reasonable description of p pho- 
toproduction cross sections on the proton, as well 
as dilepton spectra on deuterium, emerged without 
major adjustments. The key test of the spectral 
function has been provided by the dilepton ("ex- 
cess" ) spectra off iron. With average densities esti- 
mated from the decay kinematics for incoming pho- 
ton energies as used at JLAB, the rather moder- 
ate in-medium effects reported by the CLAS ex- 
periment arc fairly well reproduced. The main dif- 
ference compared to the stronger effects observed in 
heavy-ion collisions is the rather large 3-momcntum 
of the p in the CLAS data, for which the spectral 
function of Ref. [l^l predicts a significantly reduced 
broadening. Clearly, a low-momentum cut on the 
dilepton spectra would enable a critical test of the 
predicted increase in medium effects. Further con- 
straints could be obtained by analyzing absolute 
e+e~ production cross sections (e.g., the so-called 
nuclear transparency ratio) , as the in- medium spec- 
tral width of the p is directly related to its absorp- 
tion in the nuclear medium. 



*In the transport-based |29| Breit-Wigner fits in Ref. 
the extracted p-mass is consistent with the free mass. 
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^The reduction of ApjvAr to 600 MeV entails an attraction 
of ~15MeV at i)j^=(>o in the transverse p spectral function. 
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